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ABSTRACT 



The inventory control problem is examined. A FORTRAN program is 
developed, which evaluates the performance of an inventory against 
simulated demands derived from generators according to Poisson, nega- 
tive binomial, and normal probability laws. The program is tested 
using data obtained from an actual inventory control point and four 
sets of simple decision rules. The results indicate that simulation 
is a good way to test decision rules and to demonstrate the results of 
arbitrary restraints such as budgetary restrictions. 
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1 . 



Introduction. 



Not too long ago "inventory” and "wealth" were synonyms. The wealth 
of a medieval merchant was his inventory, the sum of all his goods and 
possessions. Money was one form of inventory, though a particularly use- 
ful one, since it could be exchanged readily for other forms. Interest 
was both illegal and immoral. 

In some parts of the world, even today, a marriageable daughter 
displays a sizeable portion of her father's wealth by wearing the inven- 
tory of clothing and jewelry she will take with her as a dowry. The 
jewelry is mostly in the form of coins strung together to form necklaces 
or bracelets, and their number and size are a point of family pride. 

In the southwestern part of the United States a rancher proudly 
shows his cattle to guests. It is considered poor taste, however, to 
enquire as to their number, for this form of inventory is a good index 
of his financial status. 

Most American businessmen, whether private or government, have atti- 
tudes toward inventory quite different from those illustrated above. 
Business inventory is thought to be at best a necessary evil and at 
worst "the graveyard of American business" 

Inventory takes away the freedom to use money to meet obligations 
or to engage in more profitable investment. Inventory (more precisely, 
speculation In inventory) is blamed for excesses in business cycles. In 
government, inventory represents treasury borrowings and, for this rea- 
son, is considered to be subject to interest charges. 

The purpose of inventory is to relieve supply from constant depend- 
ence on manufacturing or other commercial sources. If the customer is 
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willing to wait until the item requested can be obtained from commercial 
stocks or from production, there is no point in maintaining an inventory, 
except, of course, that buying in quantity reduces costs attributable 
to purchase order processing. But that is another story and beyond the 
scope of this paper. The purpose of inventory in the Navy is "to par- 
tially uncouple demand from commercial supply."*' The expense of inven- 
tory is the cost of service to the customer. Given that management of 
the inventory is optimal, the question of inventory size boils down to a 
definition of the degree of service desired. 

In the past the Navy has formulated budgets to maintain inventories 

expressed in terms of "months of supply". A month of supply is the 

amount of merchandise expected to be withdrawn from the inventory per 

month, based on historical experience and expressed in dollars. Our 

2 

inventories are composed largely of insurance items. Expressed in 
months of supply, they usually look outrageously large to the civilian 
leadership in the Department of Defense and in Congress, whose experi- 
ence has been largely in private business. 

In private business, the justification for inventory is profit. To 
be sure, profit in a competitive market depends, at least partly, upon 
customer good will. To the extent that good will demands, a commercial 
resale activity maintains inventory, even though that inventory may be 
uneconomical. Within this restraint, it works strenuously to weed out 

*■1 am indebted to Cdr. S. W. Blandin, S.C., U.S.N. , for this ex- 
pression as well as for many of the ideas expressed in this paper. 

•^An insurance item is a slow-moving or non-moving item which is 
retained in inventory for service or essentiality considerations. 
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slow-moving items. The penalty for low or nonexistent inventory is 
occasional loss of sales. In the Navy, the penalty may be ouch higher. 

The shortage pena lty has been defined as the cost of a com- 
pletely inoperable ship per unit of time. It is obvious, however, 
that not all shortages will render the ship completely inoperable. 
Some shortages will only partially degrade the effectiveness of 
the ship. Partial degradation is accounted for with the concept 
of essentiality ; essentiality being a relative measure of the 
seriousness of a part shortage. 

Essentiality is determined from two basic elements: (1) mission 
effect, and (2) compensability. "Mission effect" refers to the 
operational capability of the system when the part is missing, and 
"compensability" refers to the ability to make up for the shortage 
by repair, local manufacture, substitution, or cannibalization. 
Compensability is primarily a matter of time; that is, it deter- 
mines how long the mission effect loss will be felt. Judgement 
of qualified personnel is used to classify items with numerical 
values for mission effect and compensability. Essentiality is the 
product of these two numerical values. 

The shortage cost associated with a part shortage is the product 
of essentiality and the shortage penalty [ll] . 

Shortage penalties for representative ships are as follows [ll] : 



Sho rtage Pena lty 



Ship Type 


Per Year 


Per Month 


Per Day 


Guided Missile 
Heavy Cruiser 


$11,566,000 


$963,000 


$31,600 


Attack Aircraft 
Carrier (Forrestal) 


24,712,000 


2,060,000 


67,600 


Radar Picket 
Destroyer 


2,244,000 


186,500 


6,130 


Radar Picket 
Nuclear Submarine 


8,500,000 


703,000 


23,300 


Ocean Minesweeper 


690,000 


57,500 


1,890 



The above figures are based on depreciated construction and conversion 
costs and total annual maintenance and operating costs. 

Although the concept of essentiality has been understood for years, 
and although it has always been applied, on a more or less intuitive 
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basis, we have not yet been able to identify a significant number of the 
items in the Navy inventory with meaningful essentiality codes. We are, 
therefore, unable to describe the penalty cost of low inventory in a 
budget request. Essentiality is one approacn to more meaningful budget- 
ing. 

Apart from essentiality, Navy inventory managers need to justify 
budgets in terms of effectiveness, the ability to satisfy anticipated 
demands. We need to explain to Department of Defense officials and to 
Congress our need for inventory dollars in terms of performance, not 
months of supply or inventory to sales ratios. Wc need, most of all, 
some way to demonstrate to ourselves and to others the expected results 
of budget decisions. 

In the Navy as in any large organization, we often develop complex 
procedures and install them without adequate testing. Such procedures 
frequently fail because they are based on an erroneous or partial 
apprehension of reality. Or, on the other hand, the procedure may be 
correct enough and still fail because of faulty execution. Sometimes 
the wherewithal to make a procedure work is missing, whether it be funds, 
personnel, time, computer capacity, or some other essential ingredient. 

As an illustration of both kinds of failure, I offer the attempt to 
identify recurring and non-recurring demand. A demand which is non- 
recurring at a retail issue point may well be recurring when considered 
system-wide. For example, suppose a particular equipment is overhauled 
at a shipyard. As far as the shipyard is concerned, there is a one-shot 
requirement for repair parts. But the same equipment will be overhauled 
next quarter, only on a different ship at a different shipyard. Clearly 
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Che repair parts requirements are recurring, yet it is impossible to 
recognize them as such at the retail level. But even more significantly, 
the definitions of "recurring" and "non-recurring" are so unclear and 
subject to various interpretation that any conclusion drawn from such 
categorization is wishful thinking, particularly when we remember that 
it is frequently personnel at the GS-3 level who make the distinction. 

We are all familiar with the failure of some elaborate system which, 
having been installed with great fanfare, later fails decisively. We 
are fortunate when the failure is recognized, and the procedure is 
scrapped, formally what happens is that the procedure is kept because 
we are "committed" or because we cannot bring ourselves to admit the 
mistake. 

What is required is some way to evaluate decisions before we are 
totally committed. It is the purpose of this thesis to show how a 
simple simulator can help both in budgeting and in the evaluation of 
inventory control stockage rules. 
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2 . 



Inven Lory po 1 icy , 

The inventory control problem consists of two major questions: 

when to reorder and how much. In the ideal situation the demand to be 
experienced during lead time is known v as is the lead time* By lead 
time, we mean the total elapsed time between recognition of the need to 
reorder and actual receipt of the merchandise,, The reorder point is set 
at precisely the known demand corresponding to the known lead time. The 
order quantity is set by a formula which balances ordering costs against 
holding costs and calculates an "economic order quantity". At the ex- 
piration of lead time an optimal quantity resupply is delivered just as 
the last unit of the old inventory is shipped. Customer service has 
been perfect, and total expense has been minimized. 

Life in the real world of inventory management is not so predict- 
able as in the utopia described above. Demand fluctuates widely and 
unpredic tably ; lead time varies; inventory restrictions are arbitrarily 
imposed through budget restraints. Management reacts to these realities 
by striving for the best balance between service and inventory invest- 
ment * Policies are translated into decision rules which control the 
answers to the questions, when to reorder and how much. 

The U. S. Navy Ships Parts Control Center attains budget objectives 
by adjusting two arbitrary constants and by restraining economic order 
quantity and risk. Let us examine SFCC*s system more closely as an 

l 

1 Much of the advanced thinking on Navy inventory management has 
been done at the Ships Parts Control Center, Meehan icsburg, Pennsy Ivania, 
The material in this section is drawn largely from the reports written 
by the Special Assistants for Advanced Logistics Research and Develop- 
ment, of S.P.C.C. 
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illustration of scientific inventory management in the real world. 

The formula for economic order quantity used at SPCC is as follows: 



EOQ = -\1 x C 2L-1 . T 

\1 I x P 

whe re : 

Q = quarterly demand 
C = order cost 
I = holding cost rate 
P z unit price 

s a constant used for attaining management objectives 
(which, being interpreted, means living within the 
budget). 

SPCC imposes a further control by restricting the economic order quantity 
to values between one quarter's demand and five years' demand. 

At SPCC risk is defined as the acceptable probability of running 
out during lead time. Risk is calculated as follows: 



R 



Pxl x EOQ „ n 
4Q x S A 2. 



where : 




Pxl = 


holding cost 


4Q = 


annual demand 


S = 


shortage cost 


EOQ : 


economic order quantity 


A* ■ 


a multiplier 



Risk is normally restricted to values between one percent and 50 percent. 
Protection level, which is defined as one minus risk, is used in the 



7 



computation of variable safety level. 



At this stage, we must leave statistics and begin to use the tech- 
niques of probability. No matter how stable our usage rates might be - 
and usage rates in the Navy are almost never stable - any forecast of 
future usage is sure to be somewhat in error. For this reason, we base 
our reorder point on some type of projection of past demand plus a 
safety level to allow for variability. The amount of this variable 
safety level depends on the protection we desire for each line item. If 
we eliminate safety level and set reorder point at exactly the forecast 
of demand expected during lead time we run a fifty percent chance of 
running out if demands arc normally distributed about some mean or 
average value. 

The average demand is forecasted using exponential smoothing, a 
technique originated by Robert G. Brown QJ . The formula is as follows: 
Q = aD + (1-a) Q 
where : 

Q = new forecast of quarterly demand 
a * smoothing weight or constant 
D = demand experienced during past quarter 
Q = forecast for past quarter 

The beauty of exponential smoothing is that it retains the effect of 
past experience without our having to store large amounts of raw data. 
Another desirable quality is that later data are weighted more heavily 
than older data, the value of the weighting being a function of the 
smoothing constant. A further convenience lies in our ability to vary 
the smoothing constant to obtain either gradual or rapid response to 
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change, as desired. For a new item, for example, the constant can be 
set high, say at .3 to .5. As usage data arc acquired, the constant can 
be lowered, to ,1 perhaps, or even lower, so as to give less relative 
effect to last quarter's demand and more to older history. 

We are interested not only in the average demand experienced in the 
past but also in the variations from this average. Specifically we are 
interested in determining whether or not past demands fit any known 
probability distribution. If they do, we can use this distribution to 
predict the probability of our being able to satisfy forecasted demand 
with a given stock. 

By careful analysis, SPCC las found that its items fit three distri- 
butions fairly well. If the quarterly demand is .5 or less, SPCC assumes 
a Poisson distribution. If the demand is greater than .5 but less than 
or equal to 25, the negative binomial distribution is assumed. If the 
demand is greater than 25, it is assumed to be normally distributed. 

These distributions are considered further in Appendix I. 

To sec how the reorder point includes both forecasted demand and 
variable safety level let us look at a normally distributed item. 

The protection level is computed by recursively figuring the proba- 
bility that demand during lead time will be less than or equal to an 
increasing number of units of stock. When this probability builds up to 
the previously calculated protection level the process is stopped and 
the number of units corresponding to that probability is the reorder 
point less program stock and obligations. 

We see, therefore, that the inventory manager implements a budget 
restraint by scaling down the protection level through the adjustment 
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of the parameters which govern Its calculation. The effects of such a 
restraint are not felt immediately, for lead times for our technical 
items are in the order of quarters, or even years. The effect of a 
lower budget for stock replenishment will begin to manifest itself 
toward the end of the smallest lead time and will become increasingly 
evident thereafter. 




The forecasted demand protects the lower half of the demand distribution. 
The addition of variable safety level takes us up to the protection level 
previously calculated. The area to the right of the protection level is 
the risk. 
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Reorder point, then, is calculated as follows: 



RP s LQ -f VSL + P +■ 0 
where : 



L = 


lead time in quarters plus one standard deviation to 
protect against lead time variations 


Q = 


quarterly demand forecast 


P = 


program stock 


0 r 


obligations (back orders) 


VSL : 


variable safety level 
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3. Inventory Control Simulation. 

In 1962, International Business Machines Corporation completed its 
Modular Inventory Management Simulator, an elaborate computer program 
in the FORTRAN compiler language. The program is segmented so as to 
permit great flexibility. It is possible, for example, to use one of 
the forecasting routines in the program or to insert the routine used 
loci 1 y. The program includes all the functions normally associated 
with inventory management: demand forecasting, error measurement, sales, 

orders, receipts, etc. 

The program written for this thesis is an attempt, on a much smaller 
scale, to show how demand can be simulated according to probability dis- 
tributions and how such simulated demand can be used to evaluate inven- 
tory performance under various sets of decision rules. The demand gener- 
ators are subroutines, which can be used in any program where random 
variable generators are required. The program itself is reproduced in 
Appendix II. 

The first ten instructions of the main prograin accomplish the pre- 
liminary housekeeping tasks of recording the number of line items under 
consideration, setting aside storage space for the variables, zeroizing 
the distribution counters and starting the random and normal (0,1) 
number generators. 

The second segment reads in the raw data, computes the probability 
distribution parameters for each item, identifies the distribution 
applicable to each item, counts the number of items applicable to each 
distribution and calculates for each item the average quantity per 
requisition. In a normal situation these data would be continuously 
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available as one of the requirements of normal business. 



The expected demand for each item is calculated by averaging the 
first three quarters and then applying the exponential smoothing formula 
to this average and the fourth quarter demand. A high smoothing con- 
stant, 0.3, is used because the data came from a new inventory manage- 
ment center and its demands are not believed to have settled down yet. 
Similarly, for the first three quarters the squared demands are averaged. 
The expected squared demand is computed from this average and the fourth 
quarter squared demand by means of exponential smoothing with the same 
smoothing weight, 0.3. The variance is then obtained by the familiar 
formula; 



In assigning items to probability distributions it is necessary to 
test the standard deviation-to-mean ratio of items having average quar- 
terly demands of over 25. If this ratio exceeds three the item is 
identified as negative binomial distributed. This step is required 
because the probability of having a negative demand must be zero, or 
practically zero. If the mean is three tiroes the standard deviation and 
the normal distribution holds, the probability of a negative demand is 
negligible (.0035). 

The program now completes setting up the problem by recording stock 
quantities in accordance with decision rules which may be simple or very 
elaborate. These rules are in a subroutine, which can be changed at 
will. The inventory is priced, extended, and totaled. 

Finally, the program secures simulated demands from the demand gen- 
erators, extends and totals issues, and calculates effectiveness, the 
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ratio of total quantity issued on demand to total quantity demanded. 

FORTRAN is the computer language used for this thesis. More spe- 
cifically, the program was written for the FORTRAN 60 compiler for the 
Control Data Corporation 1604 computer. 

FORTRAN permits the use of symbolic machine coding interspersed 
with FORTRAN statements. Symbolic coding is used in the subroutines 
which generate pseudo random numbers and pseudo random normal (0,1) ran- 
dom variables; and its use, though highly desirable, does reduce the 
generality of the program, for the symbolic language is that of the 1604. 
Other machines have different symbolic languages, and the instructions 
would require translation before the program could be run on another 
computer. One possible solution would be to substitute generators 
written for another machine. Such generators are widely used and should 
be obtainable. 

A second limitation in these two subroutines is that both presume a 
48-bit word length. 

The Control Data 1604 general purpose digital computer at the Naval 
Postgraduate School was available for this project. This computer has 
a storage capacity of 32,768 48-bit words. 
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4. 



Demand Generators. 



The simplest of the three demand generators is that for the normal 
distribution. A subroutine was available at the Postgraduate School 
which generates pseudo random normal (0,1) variables. The number secured 
from this subroutine is converted by multiplication by the standard devi- 
ation of the item and then adding the predicted quarterly demand. The 
result is rounded off to the nearest whole number. 

For the Poisson and negative binomial items a pseudo random number 
in the range zero to one is secured from a generator, which was also 
available. This number is regarded as a cumulative probability mass. 

The probability masses for the item under consideration are calculated 
and summed thus: P(0) •+■ P(l) -hP(2) +*. . . . until the total mass equals 

or exceeds the random number. The last and highest value of the random 
variable is the simulated demand. Computationally , it is simpler to 
compute P(0) , subtract the random number, and then continue to calculate 
and add probabilities to this difference until it becomes either zero or 
positive. 

For the Poisson distribution 
P(0) = e’ q 

where : 

p(0) * the probability that the random variable is zero 
e = the natural base, 2.718 
q = the forecasted quarterly demand 
Thereafter, prolabilities are calculated recursively according to the 
formula 

p(D s 

i 
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where • 



P(X) s the probability that the randon variable is equal 
to I 

P(I-l) z the probability that the random variable is equal 
to I- 1 

q = the forecasted quarterly demand 

The iterations are terminated when I = 10, since the probability that the 

random variable will exceed 10 is negligible. 

For the negative binomial generator the parameters p,q, and r must 

be computed as follows: 

P r AM 
V 

q = 1-p 

r : AM x p 

q 

where; 

p = the probability of success for a single Bernoulli trial 
q = the probability of failure for a single Bernoulli trial 
r s the number of successes (called AR in the program) 

AM « the forecasted quarterly demand 
V s the variance in quarterly demand 

The random variable is the number of failures before the rth success. 

The probability of zero is thus: 

P(G) s p r , and thereafter 

P(I) - P( I-l) x (r + 1-1) x q 

I 
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In the iterative scheme of this program the largest value that is 
permissible for I is 16,383, this figure being the upper limit for 
indexing a "do loop". This limitation could be overcome at the expense 
of a few more instructions. 

The purpose of all the demand subroutines is to generate random 
"noise" around the demand forecasted by some manipulation of historical 
data. 
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5. Testing the Program. 

The raw data used for the test consisted of four quarters of 
demand information on 277 Items in Federal Supply Class 2815, Diesel 
Engines and Repair Parts. The data were secured from the Defense 
Construction Supply Center, Columbus, Ohio, and are typical of demand 
in repair parts classes in that high variance~to-raean ratios are preva- 
lent. The breakdown into fast, medium, and slow movers was not typical 
since 94 (34%) of the items had average demands of over 25 per quarter, 
155 (56%) had average quarterly demands of 25 or less but greater than 
or equal to 0.5, while only 28 (10%) had average quarterly demands of 
Less than 0.5. The corresponding (estimated) percentages on all cur- 
rent SPCC items are, by way of contrast, 14.8%, 26.1%, and 59.1%, 
respectively. 

It is recognized that 277 items and four quarters of usage data 
are not enough for statistical accuracy. They are enough, however, to 
show whether or not the program works. 

When the data were smoothed and checked for high variance- to-raean 
ratios, only eight of the 277 items were identified as normally dis- 
tributed, while 245 were negative binomial distributed, and 24 were 
Poisson distributed. 
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Four sets of very simple stockage rules were tested. 



Forecasted Quarterly Demand 


Unit Price 


Stock 


First Rules 
1 Poisson Items 






^ .25 


^ $100 


2 Ea 


^ .25 


$100 


3 Ea 


^ .25 


$100 


1 Ea 


.25 


$100 


2 Ea 


Negative Binomial and Normal Items 






All 


$100 


U 


All 


$100^ U/P-^SIO 


U-hS 


All 


^ $10 


U+2S 


Second Rules 

Poisson Items Same as Above 






Negative Binomial and Normal Items 






All 


Ail 


U 


Third Rules 

Poisson Items Same as Above 






Negative Binomial and Normal Items 






All 


$50 


U 


All 


^ $50 


U-hS 


Fourth Rules 
Poisson Items 






.25 


$100 


1 Ea 


Otherwise 




0 


Negative Binomial and Normal Items 






All 


^ $100 


.5U 


All 


$100-^U/P-^$50 


U 


All 


$50 — U/P^ $10 


U+S 


All 


$10 


U+2S 


Where U is the forecasted quarterly 
demand standard deviation. 


demand and S is the 


quarterly 
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For each set of scockage rules, stock levels were established. 50 
simulation runs were then made using demands from the demand generators. 
The average results are shown below. They show how effectiveness, dollar 
sales, and sales-to-inventory ratios are affected by changes in stockage 
decision rules. The evaluation of these and any other results depends 
on the goals of management. 

In the first three sets of rules, the 24 Poisson items were held at 
relatively high inventory levels. For the remaining items, various com- 
binations of unit price criteria were used to determine the si 2 e of the 
variable safety level. The second set of rules is particularly inter- 
esting in this respect because it shows the result to be expected if no 
safety level is provided. 

The fourth set of rules drastically reduced the inventory levels 
for Poisson items and applied more elaborate criteria to the remainder. 
The result was a dollar inventory of only 76 percent of the first one, 
which yielded comparable effectiveness and a better sales- to- inventory 
ratio. This result was possible because the seldom demanded Poisson 
items tended to be high priced. 



Inventory 


Sales 


First Rules 
$276,615.98 


$143,501.35 


Second Rules 
$182,009.56 


$118,389.19 


Third Rules 
$239,530.38 


$136,101.93 


Fourth Rules 
$211,419.04 


$118,329,89 



Sales/Inventory 


Effectiveness 


e 

00 


.952 


,651 


.680 


.569 


.878 
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.561 



.947 



6. Concilia loots . 



In significant positions * notably at the top - the Department of 
Defense is staffed with civilian personnel who understand mathematical 
models, gaming, and probability theory. These people rightly insist on 
discussing budget problems in such terms. It behooves us in the military 
to learn to speak their language. Inventory control simulation with 
probabilistic demand generation is one way to do so. 

To be useful to an inventory control point, a simulator needs to be 
more elaborate than the one written for this thesis. It should be con- 
sistent with the other routines of the center, so that it would be pos- 
sible, at any time, to leave off normal processing and project as far 
into the future as is desired, using simulated demand and, if desired, 
revised decision rules. It should take into consideration both antici- 
pated receipts (dues) and back orders (obligations). It should establish 
dues and obligations and should record receipts and the release of obli- 
gations. 

The IBM Modulator Inventory Simulator does all these things. To- 
gether with probabilistic demand generators, it could form the basis of 
a local program. 

The Navy is always faced with the task, not only of understanding 
its problems, but of communicating them to the civilian leadership and 
to Congress. Inventory control simulation can help us do both. It is, 
first of all, a good way to look at our own business and at proposed 
changes. It is, further, a way to present budgets, showing clearly the 
anticipated results of various levels of funding. It deserves serious 
consideration as a management tool. 
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APPENDIX I 



PROBABILITY DISTRIBUTIONS 

The normal probability Law is too well known to require treatment 
in this paper. It is thoroughly discussed in almost any probability 
text. 

Less familiar is the Poisson distribution, also called the law of 

rare events. It is defined as: 

The number of occurrences of events of a specified type in a 
period of time of length 1 when events of this type are occurring 
randomly at a mean rate per unit time [fj . 

The Poisson probability mass function is 
- ^ "2 x 

e x = 0,1,.. . 

x ; 

Its mean and variance are both equal to /L • 

Even less familiar is the negative binomial or Pascal distribution. 

This distribution is defined as: 

The number of failures encountered in a sequence of indeDendent 
Bernoulli trials (with probability p of success at each trial) 
before the rth success [V] . 

The negative binomial mass function is 




Its mean is ^ while the variance is — . The negative binomial 

P p 2 

fits the demand distribution of many items of supply. One of its desir- 
able properties is that, while the probability of the random variable 
being equal to zero may be quite high, it is impossible for it to be 
negative. This is a characteristic which makes the negative binomial law 
fit demand experience of items having relatively high variance to mean 
ratios, where the normal distribution would be inaspropriate . 
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PROGRAM DATA! 

N = ?77 

DIMENSION GUSAGE! 300 ) , VAR ( 300 ) , REQAVE ( 300) , 

1 ISTKNRI 300) , ATI ME ( 300 J , S I GMA ( 300 ) , 

2 IND I C ( 300 ) .PR ICEI/300) , STOCK! 300 ) , 

3DEMAN0! 300) , S TK EX T ( 300 ) , USEXT! 300) 

COMMON QUSAGE, STOCK, INDIC, PRICt, SIGMA 

NUNIF=1695 

I P I N0M = 0 

NEGD I N=0 

IP0IS=0 

CALL RANDOM ( 0 , R ) 

SETUP PECULIAR TO TEST DATA 
00 15 1=1, N 

1 ORE AD 2 I STKNR! I ), PRICE! I ) , IREC2, IQREC2. 
1NREC2,IQNREC2, l RE C3 , 1 QREC 3 , NREC3, IQNREC3, 
2IREC4, IQRECU.NRECU, IQNREC4 , I REC 1 , IQREC1 , 

3NREC1 , IQNREC1 

2 FORMAT! 1 7 , 1 2X , F 7. 2 , 5X , 8 ( 1 2 , 1 4 ) > 

NR 1 = I REC 1 + NREC1 

‘NR2 '= IRLC2 +' NREC2 
NR 3 = IREC3 + NREC3 
NR4 = I R tC4 + NREC4 
QTYl = IGREC1 + IGNREC 1 
QTY2= IQREC2+I GNREC2 
QTY3= I QREC 3+ I UNRE C3 
QTY4=IQREC4+ICNREC4 
NR = NR 1 + NR2 + NR3 + NR 4 
QTY=QTY1+QTY2+QTY3+QTY4 
QUSAGE! I ) = (QTY2+QTY3+QTY4)/3. 

AVESQ = (QI Y2*Q1Y2+QTY3*QTY3+QTY4*QTY4) /3. 

ALPHA= . 3 

QUSAGE! 1 )=ALPHA*3TY1+! 1 . -ALPHA ) *QUS AG E ( I ) 
AVESQ=ALPHA*OTY 1* QT Y 1 + ( 1 . - ALPHA » * A VES Q 
CALCULATION OF VARIANCE AND STANDARD DEVIATION 
VAR! I ) = AVESQ-QUSAGE( I )*QUSAGE ( I ) 

S I GMA ( I ) = SQRTF ( VAR ( I ) ) 

IDENTIFICATION OF DEMAND DISTRIBUTION 
IF! SIGMA! I ) -QUSAGE ( I )/3. ) 8,8,4 
8 IF (QUSAGE! I ) -25. ) 4,4,3 

3 I ND I C ( I ) = 1 



7 

15 
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IBINOM= IBINOMtl 
GO TO 7 

I F ( QUSAGE ( I )- . 5 ) 6,6,5 

INDIC ( I )=2 

NEGBIN=NEG8IN+1 

GO TO 7 

INDIC! I ) =3 

IPOIS=IPOIS+1 

ANR = NR 

REQAVE! I ) =QTY/AN3 
STKSUM=0. 

CALL BUY 1 ( N ) 

EXTEND INVENTORY AND TOTAL 
D04 1 1=1. N 

STKEXT! I 1 = PRICE ( I )*STOCKII) 
STKSUM=STKSUM+STKEXT{ I ) 
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c get demand from demand generators and extend 

USET0T=0. 

DEMEXT=0. 

ANUM=0. 

OENOM = 0 » 

004 £L I = 1 , N 

I F ( INC I C ( I ) -2 ) 42,43,44 

42 CALL NORMAL ( NUNI F , QUSAG E ( I) , S I GMA ( I ) , AN ) 

GO TO 45 

43 CALL NEG8INO(QUSAGE(I ),VAR(I),AN) 

GO TO 45 i 

44 CALL PO I SSON( QUSA GE ( I ) , AN) 

45 OEMANC ( 11= AN 

DEMEXT = DEMEXT+DE'1AND( I )• PRICE! I ) 

C CALCULATE EFFECTIVENESS 

l F ( STOCK ( I )-DLMAvJO ( I ) ) 49,46,46 
49 ANUM=ANUM+STOCK(I ) 

US EXT ( I ) =STOCK( I) * PRICE! I ) 

GO TO 47 

46 ANUM= ANUM + DEMAND! I ) 

US EXT ( I ) =DEMAND ( I >* PR ICE Cl) 

47 USETOT=USETOT+USEXT( I ) 

48 DENOM=DENOM+DEMA'JD ( I ) 

EFFECT=ANUM/DENOM 

52 PRINT 51 STKSUM,JSETOT, EFFECT 
51 FORMAT! IX, 3IF15.3 ) ) 

END 



> 
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C POISSCN CEMANO GENERATOR 

SUBROUTINE PO I SS3N ( Q , AN) 

1 CALL RANDOM! 1,R) 

2 PROB=2.7 18**(-G) 

3 SUM=PROB-R ! 

4 I F ( SUM ) 7,5,5 

5 AN=0 . 

6 GO TO 13 

7 DO 12 N= 1 , 10 
B AN=N 

9 PROH=PROR»Q/AN 

10 SUM = SUM+ PROB 

11 IF! SUM) 12, 13, 13 

12 CONTINUE 

13 END 
C 

C NEGATIVE BIN0MIA W DEMAND GENERATOR 

SUBROUTINE NEGBINO l AM, V, AN) 

1 P = AM / V 

C P RESTRICTED TO VALUES BETWEEN .000001 AND .999999 

l IF(P-. 000001 ) 202,2,200 

200 IF(P-. 999999) 2,2,261 

201 P=. 999999 
GO TO 2 • 

202\p= . 00000 1 

2 G^= 1 .-P 

3 AR=AM*P/C+.5 

5 I AR= AR 

6 AR= I AR 

7 IF! AR ) 8,B, 10 

8 AR= 1. 

10 CALL RANDOM ( 1 , R) 

11 PROB = P* » AR 

12 SUM = PROB - R 

13 I F l SUM ) 16,14,14 

14 AN=0 • 

15 GO TO 22 

16 DO 21 1=1 , 16383 

17 AN= I 

18 PROB=! AR + AN-1 . ) MN*PROB*Q 

19 SUM = SUM + PROB 

20 I F ! SUM ) 21,22,22 

21 CONTINUE 

22 END 
C 

C NORMAL DEMAND GENERATOR 

SUBROUTINE NORMAL ( NUNIF, AM, S, AN) 

1 CALL RNOEV (NUNIF, DEV) 

3 AN=S*DEV+AM+. 5 

N=AN v 

AN = N 

4 END 



30 



on 



; NORMAL (C,l) RANDOM VARIABLE GENERATOR 

SUPROUTINE RNUEVI NUN IF, DEV) 

CONI NZ=777777 7777 777777B) . 

RS V ( ERA= 1 I . 

ENA(O) STA('ERA) 

SILKS) F.NI1 (15) 

8RDV LDA(NUNIF) MU If 1220703125) 

SCL I 7600000000000 OOOB ) STA( NUNIF ) 

RAO ( ERA ) I J PI { 8RDV ) 

SUBI 2000000000000000B ) 

AJP2IL+1) LAC(ERA) 

ENAl(O) LLSI36) 

SCAK 2057B) LRSI47) 

SSK(ERA) SCM(NZ) 

SCM(NZ) FMUI 1724673317272054B) 

8 STA(I)EV) CN II (N) . 

RETURN 
END 



STAI ERA) 



RANDOM NUMBER GEMERATOR 
‘SUBROUTINE RANDOM (K,R) 

1 I F ( K ) 5,2,5 

B 2 E= 4000 0000 0000 0000 

B 3 E2=20C0 COOO 0000 0000 

B 4 X = 1777 7777 7777 7777 

5 SIU6IL+7) ,ENI6(2) , ENQ(N) , LDA ( X ) , L L S ( 1 91 , 
SCL(E),ADD(X),SC- (E), IJP6 ( L- 1 ) , STA ( X) . 
ARS ( 1 1 ), ADDIE2) ,F ADI E2 ) , STA ( R ) , ENI 6 IN ) , 
ENI IN) 

6 END 
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C HIRST SET OF STOCKAGE RULES 

subroltine euYim 

10DI MENS ICN GUSAGE! 300 ) , STOCK ( 300) , 

1 INDIC130C) , PRICE! 3CC) , SIGMA (300 
COMMON CUSAGE,ST3CK, I NO I C, PR ICE, SIGMA 
0020 1=1, N 

C IDENTIFICATION OF POISSON DISTRIBUTED ITEMS 

IF ( INCIC ( I )-2) 13*10,2 
C POISSCN RULES 

2 I F ( GUSAGE ( I )- . 25) 6,6,3 

3 I F ( PR ICE ( I ) — 1 00 • I 5,5,4 

4 STOCK ( I ) =2. 

GO TO 20 

5 STOCK ( I ) =3. 

GO TO 20 

6 IFCPR ICE< I )-1CC.) 8,8,7 

7 STOCK ( I ) = 1 • 

GO TO 20 

8 STOCK ( I ) =2. 

GO TO 20 

C ‘NORMAL AND BINOMIAL RULES 

10 IF C PR ICE ( I ) - 1 00 . ) 12,12,11 

11 ISTCCK = G USAGE! I K • 5 
GO TO 19 

12 IFCPR ICEC I 1-10. > 14,14,13 

13 ISTCCK = 0 USAGE! I H- SIGMA! I ) +.5 
GO TO 19 

14 I STOCK = GUSAGE ( I )«• 2 . *S IGM A ( I ) 4 . 5 

C ROUNC OFF TO NEAREST WHOLE NUMBER 

19 STOCK ( I ) = I S TOCK 

20 CONTINUE 
25 END 

S E CON C SET OF STOCKAGE RULES 
SUBROUTINE BUYl ('4 ) 

1 0D I MENS I ON QUSAGE! 3C0) , STOCK (300), 

1 INDIC (30C ) ,PR ICE! 300 , SIGMA (300) 

COMMON QUSAGE, STUCK, I NO I C , PR I C E , S I GMA 
0020 1=1, N 

C IDENTIFICATION OF POISSON DISTRIBUTED ITEMS 

IF! INCIC! I)-2 ) 10,10,2 

C POISSCN RULES 

2 I F ( QU S AGE ! I )- . 25) 6,6,3 

3 IF! PRICE! D-100.) 5,5,4 

4 STOCK ( I ) =2 . • 

GO TO 20 

5 STOCK! I ) =3. 

GO TO 20 

6 I F ( PR ICE! I ) - 1 00 . » 8,8,7 

7 STOCK { I ) = 1 . 

GO TO 20 

8 STOCK ( I )=2* 

GO TO 20 

C NORMAL AND BINOMIAL RULES 

10 ISTCCK = QUSAGE ( I )♦ .5 
C ROUND OFF TO NEAREST WHOLE NUMBER 

17 STOCK ( I ) = ISTOCK 
20 CONTINUE 
25 END 
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THIRD SET OF STORAGE RULES 
SUBROUTINE BUY1 ( N ) 

1 OC I MENS I CN GUSAGE! 3CC ) , STOCK (3CC) , 

1 IND l C ( 30C ) , PR. I C E( 3C0) , SIGMA ( 300) 

COMMON GUSAGC.ST3CK, 1NDIC, PR ICE, SIGMA 
C020 1=1, N 

IDENTIFICATION OF POISSON DISTRIBUTED ITEMS 
IF! INCIC! I ) -2 ) 10,10,2 
POISSCN RULES 

2 IF ( GUSAGE ( I )- .25) 6,6,3 

3 I F ( PR ICE ( I )-100.) 5,5,4 

4 STOCK! I ) =2 . 

GO TO 20 

5 STOCK I I ) = 3. 

GO TO 20 

6 IF! PRICE! I ) - 1 00 . ) 8,8,7 

7 STOCK! I ) = 1 . 

GO TO 20 

8 STOCK ( I ) =2. 

GO TO 20 

4 NORMAL AND BINOMIAL RULES 

10 IF! PRICE! I )-50. ) 11,11.12 

11 ISTOCK = GUSAGE( I I+SIGMA! 1 ) + .5 
GO TO 17 

12 ISTCCK = QLSAGE( I )+ .5 

ROUND OFF TO NEAREST WHOLE NUMBER 
17 STOCK ( I ) = ISTOCK 
20 CONTINUE 
25 END 

FOURTH SET OF STDCKAGE RULES 
SUBROUTINE BUY1 (N ) 

100 1 MEN SI ON GUSAGE! 300) , STOCK (300) , 

1 IN0IC(30C) .PRICE! 3CC) .SIGMA (300) 

COMMON GUSAGE, STOCK, INDIC, PRICE, SIGMA 
D020 1 = 1 , N 

IDENTIFICATION OF POISSON DISTRIBUTED ITEMS 

IF! INCIC! I )-2) 13,10,2 

POISSCN RULES 

IF! PR ICE! I J-100.) 2,6,6 

2 IF ( GUSAGE ! I J-.25) 6,6,3 

3 STOCK! I )=1 . 

GO TO 20 

6 STOCK! I )=0. 

GO TO 20 

NORMAL AND BINOMIAL RULES 

10 IF(PR ICE! I ) - 1 00 . ) 12,12,11 

11 I S TOCK= . 5*QUS AGE! I ) + .5 
GO TO 19 

12 IF! PR ICE! I >-5C. ) 14,14,13 

13 ISTOCK = QUSAGE( I )+ .5 
GO TO 19 

14 IFIPRICE! I )— 10« ) 16,16,15 

15 ISTOCK=OUSAGE( I )+ SIGMA! I )+.5 
GO TO 19 

16 ISTOCK=QUSAGE ( I H 2 .*S I GM A ( I ) + . 5 

19 STOCK! I )=ISTOCK 

20 CONTINUE 
25 END 
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